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DIFFERENTIAL  CROSS  SECTION  AND  RELATED  INTEGRALS 
FOR  THE  MOLltRE  POTENTIAL 


INTRODUCTION 

At  the  core  of  all  radiation-damage  and  ion-range  calculations  for  heavy-ion  beams 
incident  on  bulk  materials  is  the  choice  of  potential  that  represents  the  interaction  between  the 
incoming  ions  and  the  lattice  atoms.  These  potentials  range  from  ones  specifically  tailored  to  an 
atom-ion  pair  to  less  accurate  forms  that  can  be  used,  with  simply  defined  parameters,  for  any 
interacting  pair  (1],  Most  of  the  forms  in  this  latter  category  consist  of  some  approximation  to 
the  Firsov  form  of  the  two-body  Thomas-Fermi  interaction. 

The  Thomas-Fermi  potential,  for  an  isolated  atom  of  charge  Z 2e,  is  usually  written 

y(r)  - yy  Xrf'V*)- 

where  XrOO  is  the  screening  factor  for  the  Coulomb  potential  and  where  the  screening  radius 


a“2  irl  —I  ZiUi  “ 0 8853  aoZiUi- 

32  J me 1 

The  function  xtM  is  available  in  tabular  form  [la].  Firsov  was  able  to  justify  the  adaptation 
of  the  Thomas-Fermi  potential  as  a two-body  interaction;  specifically,  if  the  quantities  Z,  and 
Z2  are  the  atomic  numbers  of  the  incoming  ion  and  lattice  atom  respectively,  we  write 


K(r)  - 


Z,Z2eJ 


Xr(r/a). 


where  we  adopt  the  screening  radius  of  Lindhard,  Nielsen,  and  ScharfT  (LNS)  [2]  given  by 

a - 0.8853  a0Z~xl\ 


z - jz,3'3  + zf/3)3 


In  addition  to  the  Thomas-Fermi  screening  function  being  available  in  tabular  form,  it  has 
been  approximated  by  a large  variety  of  analytical  forms.  Of  these  forms,  we  are  particularly 
interested  in  the  Molifcre  form  [lb],  given  by 

XmM  - 0.35  e~°3x  + 0.55  e~l  i*  + 0.10 

where  x — r/a.  The  Molilre  screening  factor  falls  off  exponentially  with  large  separations, 
whereas  the  Thomas-Fermi  screening  factor  falls  off  as  x-3.  It  has  been  shown  however  that 
the  Thomas-Fermi  interaction  falls  off  too  slowly  and  that  the  Molifcre  potential  is  a more  real- 
istic interaction  for  large  separations  [lb, 31.  The  Moli&re  interaction  is,  for  example,  used  by 
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Robinson  and  Torrens  in  their  computer  simulation  studies  of  radiation  damage  (4].  In  Fig.  1 
we  show  the  ratio  of  the  Molifere  and  Thomas-Fermi  screening  factors  as  a function  of  separa- 
tion. The  agreement  is  quite  good  out  to  the  region  where  the  exponential  decay  of  the 
Molifcre  screening  factor  is  dominant.  In  this  region,  as  we  said,  the  Molfere  interaction  is 
more  realistic. 


O 2 0 3 0.5  07  1 2 3 5 7 

REDUCED  SEPARATION  x 


Fig.  1 — Ratio  of  Molifcre  and  Thomas-Fermi  screening 
factors  as  a function  of  dimensionless  separation 


Although  computer  simulation  calculations  use  an  interaction  directly,  Boltzmann  [5]  and 
Lindhard  [2,6]  transport  calculations  use  the  scattering  cross  section  associated  with  the  poten- 
tial. With  the  definitions 
p — impact  parameter, 

£ — incident  ion  energy, 
and 

A 2 

Ei  — — — £ — center  of  mass  energy, 

Ai  + A2 

where  A\  and  A2  are  the  masses  in  atomic  units  of  the  incoming  ion  and  lattice  atom  respec- 
tively, the  scattering  angle  is  [7] 

* “ W ~ 2Pln  r2[l  - V{r)/Ei-pVr*\'"'  0) 

where  rmin  is  the  largest  zero  of  the  radical  in  the  integrand.  We  also  make  use  of  another  vari- 
able, the  energy  transferred  in  a collision, 

T - Tm  sin2  »/2, 

where 


Tm  - yE 


* ^1^2 
(^i  + A j)2 


is  the  maximum  kinetically  allowed  energy  transfer.  Corresponding  to  this  maximum  energy 
transfer  is  the  minimum  energy  that  can  be  carried  away  by  the  incident  ion,  given  by 
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f min  “ E/P 


^1  + ^2 


Much  of  the  work  of  Lindhard  and  his  coworkers  is  couched  in  terms  of  the  dimensionless  vari- 
ables 


and 


« - E/El, 

Z \Z 2?^  A \ •¥  A 2 


i - f2  T/Tm  - «2  sin2  0/2. 


In  terms  of  these  various  definitions,  the  differential  cross  section  is 
da  — — 2 irpdp  — - 2 irp(t)  dt. 


LNS  use  the  notation 


so  that 


da 


f (r'^j  - ± t3"pU) 


dp(t ) 
dt 


(2) 


Out  of  a desire  to  create  a simple,  universal  cross  section,  LNS  now  make  two  approxima- 
tions to  obtain  / (rl/2):  they  replace  Eq.  (1)  by  the  momentum  approximation  [8]  to  the 
scattering  angle,  so  that 

0--p/E]  J dr  i (I  -pVr2)-'11.  (3) 

' r ar 

and  they  make  the  substitution 

t “ j^J2  ~ **  s^2  ®/2- 


By  combining  these  approximations,  we  find  that  we  can  write 


tui  “ - y ( p/a)  J dx  x 1 


,.i£^ 


-1/2 


d_ 

dx 


7 


(4) 


Equation  (4)  provides  a functional  relationship  between  the  impact  parameter  and  the  reduced 
energy  transfer  that  can  be  used  to  solve  Eq.  (2)  for  /(rl/2).  The  advantage  of  the  LNS 
method  is  that  the  differential  cross  section  depends  on  only  one  variable;  in  other  words,  the 
variable  t in  Eq.  (4)  would  ordinarily  depend  on  both  p/a  and  «,  but  in  the  LNS  approximation 
it  depends  only  on  p/a. 


Before  we  continue,  one  other  point  should  be  made.  With  any  infinite  range  potential, 
the  total  cross  section  diverges  as  the  energy  transfer  approaches  zero.  To  bypass  this  difficulty. 
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we  cut  off  the  allowed  energy  transfer  at  a minimum  value  T j.  The  cutoff  value  will  be  of  the 
order  of  20  eV,  dependent  on  the  target  material,  and  can  be  thought  of  as  related  to  the 
minimum  energy  required  to  displace  an  atom  from  its  lattice  site.  Collisions  in  which  an 
energy  of  less  than  7",  would  be  transferred  to  a lattice  atom  are  not  allowed.  Roughly  the  same 
approach  is  taken  in  computer  simulation  calculations,  where  there  is  a maximum  allowed 
impact  parameter. 

Given  these  considerations,  the  total  (macroscopic)  cross  section  is 

c 

No-(E)  - N J*  da  - Nira2  J*  dr\  t>~2  f{r\), 

^0 


where  N is  the  number  density  of  lattice  atoms  and 

f e r. 


*!o  ' 


yEl 


1/2 


Following  LNS,  we  can  write 


N S(E)  - N it  a 


di\  /( r|) 


for  the  stopping  cross  section  for  elastic  collisions  and 


N W(E)  - Nira2 


yEl 


c 

J*  d t)  I)2  /(■»)) 


for  the  square  fluctuation  in  energy  loss.  The  quantities  /,  S,  and  IF  are  provided  by  LNS  in 
tabular  form  121.  Manning  has  generated  a more  complete  table  [91.  Winterbon,  Sigmund,  and 
Sanders  (WSS)  [6]  have  created  a more  convenient  fit  to  the  LNS  /(■»>)  in  the  form 

[1-3/2 

1 + (2Xt>4/,)2/3]  , X-  1.309. 

Relatively  convenient  expressions  can  be  found  for  a(E),  S(£),  and  W(E)  [10,111. 

It  is  our  intention  in  the  rest  of  this  report  to  find  the  /(r>)  corresponding  to  the  Moliere 
potential  and  to  develop  convenient  forms  for  /,  <r,  S.  and  W. 


EXACT  RESULTS 

The  integral  that  arises  in  Eq.  (3)  with  the  Molifere  potential  is  evaluated  by  Lehmann  and 
Leibfried  [8],  with  the  result 

3 

0 - a'1  £6,X,*,(X,p/fl). 

/-i 

where  A,  is  the  modified  Bessel  function, 

X,  -0.3.  1.2.  6.  / - 1.  2.  3. 

and 

b,  - 0.35.  0.55.  0.10.  / - 1.  2.  3. 
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Consequently  we  have 


4' 


V - l'n  - 'I  />,X,AT  ,(X,/»/o). 
1 /-I 


In  terms  of  the  variable  tj,  Eq.  (2)  becomes 

« 2 17 2p/a 

and  we  find 

A — 1 

IV,2 


djL 


JjH 1 3 


dp  la 


/-1 


dp /a 
K |(X,p/o) 


\,p/a 


+ K0(X.,pla) 


(5) 


(6) 


(7) 


By  using  Eqs.  (5),  (6),  and  (7),  we  can  form  a table  expressing  the  relationships  between 
/( v)>  V and  pi  a. 

For  convenience  let  us  make  the  definitions 


*aM(£)-NwaJ[e„(T,o)  -£„(«)], 

00 

- / dri  tj_j /«(*?). 

X 

N S„(E)  -Nna*1^  [$„(«)  - SM(W|. 

X 

SM(x)  - f dri  fM( T)), 

0 

N WU(E)  -Niro1  |^|  [ o)j, 

X 

®A/(*)  ~ f dr)  tj2/w(tj). 


The  first  function  is  trivial  to  evaluate,  so  that 
The  second  quantity  is 


(8) 


(9) 


(10) 


and 


00  f 3 2 

“ \ f dq  q ]£  b,ktKt{k,q)  . 

i p(x)/a  [,-l 


This  integral  can  be  evaluated  exactly  [12],  yielding 


- j I *,2  i,: 


*«?(*,)  + - Ar0(s,)X  ,(j,)  - *,’(*,) 


T £ TT^TF  KM  - s,AT0(i,)Arl(s/)|. 

* ijm  1 V * * 


i+! 


where  sKj  - k.,jp/a.  The  quantity  fflw(jr)  apparently  cannot  be  evaluated  in  closed  form. 
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Figure  1 indicates  the  deviations  of  the  Molifere  potential  from  the  Thomas-Fermi  poten- 
tial. Specifically  we  note  the  series  of  wiggles  in  the  ratio.  These  are  no  doubt  due  to  the 
fitting  of  the  Moli&re  form  to  the  Thomas-Fermi  potential;  they  have  no  physical  significance. 
In  Fig.  2 we  reproduce  that  curve,  with  the  addition  of  a plot  of  the  ratio  of  the  Molifere  and 
Thomas-Fermi  cross  sections,  drawn  as  a function  of  p/a.  We  see  that  the  wiggles  in  the 
Motifere  potential,  relative  to  the  Thomas-Fermi  potential,  are  reproduced,  and  indeed 
magnified,  in  the  corresponding  cross  sections.  In  part  this  magnification  is  probably  due  to  the 
use  of  the  momentum  approximation  to  calculate  the  scattering  angle.  Figure  3 provides  plots 
of  the  Moliere  and  Thomas-Fermi  f(i j)  functions;  also  shown  is  the  WSS  fit  to  the  Thomas- 
Fermi  f(y)  function. 


Fig.  2 — Ratio  of  the  Moliere  and  Thomas-Fermi 
screening  functions  as  a function  of  r/fl,  and  ratio  of 
their  /ft •/*)  functions,  as  a function  of  p/a 


SIMPLE  FIT  FOR  /*( tj) 

It  is  time  consuming  to  evaluate  the  exact  expression  (Eqs.  (5),  (6),  and  (7))  for  /*/(*}); 
in  actual  use  it  is  convenient  to  have  some  simple  form  to  represent  the  cross  section.  We 
have  found  that 


/(t|)  - 

“ fjiv) 


a0y  In  t)  + a|T|  + a2r|2  + ajy3,y  < y* 


01  + V2 

7,  T)  > T) 

02  + 0JTl  + V2 


0.06. 


fits  the  exact  to  better  than  6%  for  all  values  of  y.  The  values  of  the  parameters  in  these 

equations  are 


a,  - -20.45,  -71,  422.097,  -1429.70,  / - 0,  1,  2.  3, 
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REDUCED  ENERGY  TRANSFER  i)*tl/z 

Comparison  of  the  Moliere,  Thomas-Fermi, 
and  WSS  scattering  functions 


and 

0,  - 0.007.  0.0387,  0.826).  / - 1.  2,  3. 

In  particular,  a2  and  a2  were  chosen  so  that  the  fits  in  the  two  regions  had  equal  values  and 
derivatives  at  ij  - r)  *. 

The  various  integrals  of  /( 17),  Eqs.  (8),  (9),  and  (10),  can  be  evaluated  in  closed  form. 
For  t|  < t;  * we  have 

C 1(17)  — 148.298  — y a0  In2 17  + a.  In  17  + a2T)  + y ajT72|, 

Si(-»|)  - y a0T}2  |ln  17  - yj  + y a,i72+  y a2T)3  + y a3l74. 

and 

©1(9)  - y i74a0|ln  v - y + y <*\V*+  y aw5  + y 
and  for  tj  > sj  * we  have 

“ {[(0203  ~ 20,(0/  - 2 02)]/>2(i7)  + 4/3,^2/t, 

+ (02  - 2/Ji03)[p1(t,)  - 2 In  tj)]J/(40/), 

S2(t>)  - 0.059298  + y [(4/8,  - 0j)/’2(t»)  + />1(tj)J. 
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and 


©2(17)  “ 0.079916  + j ji,2  + (4/9,  - 2/9^tj  + (/932  - /92  - 2/91/93)P1(t,) 
+ [(2/9,Of  - 2/3 2)  - /93(/932  - 3/92)Jp2(tj) j. 


where  we  write 


and 


P ,(•»})  " In  (/3j  + t>/93  + if , 

Pi(v)  - Pq'  In  [(/93  + 2if  - />o)/(/33  + 2tj  + P<>>]- 


Fig.  4 — Comparison  of  present  work  with 
exact  Molifcre  results:  ratio  of  scattering  ker- 
nels as  a function  of  reduced  energy  transfer 


Figure  4 presets  a comparison  of  our  fit  (/,  and  /2)  with  the  exact  expression  /m(tj). 
We  see  that  the  error  in  /(tj)  is  never  more  than  6%.  Consequently  the  errors  in  C,  $,  and  © 
will  be  less  than  6%.  Because  we  are  now  prescribing  an  /(tj),  we  are  interested  in  what 
screening  factor  corresponds  to  this  new  /(tj).  Within  the  spirit  of  the  LNS  approximations, 
embodied  in  Eq.  (4),  we  can  write  [lc,  13,  14] 


X(x)  - — f dq  (?2  - x2)  r\(q) 

n X ' ' 


where  x and  t>,  are  related  by 


- J dtf  7,  2/( tj). 


In  Fig.  S we  compare  as  a function  of  x the  *(x)  obtained  in  this  manner  from  our  new  /(tj) 
with  the  exact  form  x*(x).  We  see  that  *(x)  differs  from  *M(x)  by  less  than  6%.  We  also 
show  the  error  in  /( tj)  again,  this  time  as  a function  of  q — p(i\)/a. 
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x = r/o  AND  q = p(^)/o 


Fig-  5 — Comparison  of  present  work  with  exact  Molifcre  results: 
ratios  of  scattering  kernels  and  screening  functions 


CONCLUSIONS 

We  have  created  a simple  differential  cross  section  that  reproduces  the  cross  section 
derived  from  the  Molifcre  potential  by  the  LNS  method.  The  expressions  for  this  /( tj)  and  the 
related  integrals  for  total  cross  section,  stopping  power,  and  fluctuation  in  energy  loss  involve 
only  simple  powers  and  logarithms.  We  suggest  the  use  of  this  version  of  the  Molifere  cross 
section  as  a universal  cross  section,  to  replace  the  LNS  and  WSS  forms. 
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